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Abstract 

The B® — > J/ipKs decay has recently been observed by the CDF collaboration 
and will be of interest for the LHCb experiment. This channel will offer a new 
tool to extract the angle 7 of the unitarity triangle and to control doubly Cabibbo- 
suppressed penguin corrections to the determination of sin 2/3 from the well-known 
B® — > J/ipKg mode with the help of the {/-spin symmetry of strong interactions. 
While any competitive determination of 7 is interesting, the latter aspect is partic- 
ularly relevant as LHCb will enter a territory of precision which makes the control 
of doubly Cabibbo-suppressed Standard-Model corrections mandatory. Using the 
data from CDF and the e + e~ B factories as a guideline, we explore the sensitivity 
for 7 and the penguin parameters and point out that the B®-B® mixing phase cj) s , 
which is only about —2° in the Standard Model but may be enhanced through new 
physics, is a key parameter for these analyses. We find that the mixing-induced 
CP violation S(Bg — > J/ipK§) shows an interesting correlation with sin</> s , which 
serves as a target region for the first measurement of this observable at LHCb. 
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1 Introduction 



As the LHCb experiment at CERN's Large Hadron Collider (LHC) has just started 
to take first physics data, we are entering a new territory in precision flavour physics. 
A particularly exciting feature of the LHCb research programme is the exploration of 
decays of B s mesons. Among the various B s decays with a promising physics potential 
is the B° s — » J/ipKs mode, which has recently been observed by the CDF collaboration 
at the Tevatron [I]. 

InRef. |2], it was pointed out that the CP-violating asymmetries of this channel offer 
a probe for extracting the angle 7 of the unitarity triangle (UT). Moreover, it was noted 
that also penguin topologies can be determined, which limit the theoretical precision of 
the measurement of sin 2(3, where (3 is another UT angle, by means of the mixing-induced 
CP violation in the "golden" decay B® — > J/ipK^. In this strategy, the [/-spin flavour 
symmetry of strong interactions is used, which is a subgroup of S77(3)f relating down 
and strange quarks to each other, and allows us to relate the hadronic parameters of the 
B® — > J/ipKs and B® — > J/ipKs modes to one another. In Fig. [IJ we show the decay 
topologies of the B® — > J/ipK s channel; interchanging all down and strange quarks, we 
obtain the B® — > J/ipK s topologies. 
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Figure 1: Decay topologies contributing to the B° s — > J/ipK s channel: "tree" (left) and 
"penguin" (right) diagrams. 

In view of the start of LHCb, we would like to have a closer look at the prospects 
for analysing the B® — > J/ipKs channel at this experiment. Here we will give a special 
emphasis to the extraction of the hadronic penguin parameters and the control of the 
corresponding uncertainty in the measurement of sin 2/3 from B® — > J/ipKo,. In Ref. |3), 
using data on the B® — > J/ipir channel (see also [i]), it was pointed out that such 
effects tend to lower a tension in the fit of the UT between measurements of its side 
Rb oc |V^b/V^(,| and sin 2/3. Following these results, we will explore the prospects for such 
an analysis by means of the B® — > J/ipKs channel as an analysis of B® — > J/ipir is 
considered to be more challenging at LHCb. 

Due to the low branching fraction of the B° s — > J / ipKs mode, LHCb might be able to 
get only very crude first information on the corresponding CP asymmetries by the end 
of 2011, even though a signal should be clearly visible by then. In order to explore the 
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full potential of this decay, we will therefore assume two milestones for our feasibility 
study: 

• an integrated luminosity of 6 fb _1 , which the LHCb experiment could realistically 
collect by the end of the second data taking period of 2014-15 (with the LHC 
having reached the design centre-of-mass energy ^Js = 14 TeV); 

• an integrated luminosity of 100 fb -1 , which is the current target for a subsequent 
LHCb upgrade. 

For the measurement of sin 2/3 through B® — > J/tpK s , (sin 2/3) j/,^, LHCb expects 
precisions of 0.022 with 2 fb" 1 |], and 0.003-0.010 for 100 fb" 1 , depending on to be 
determined systematic errors [6]. In view of these impressive numbers, we have to take 
penguin effects into account in order to match the experimental with the theoretical 
precision. Such an analysis may eventually also allow us to pin down new-physics (NP) 
effects in B®-B® mixing; the current data still leave space for NP contributions as large 
as about 50% of the SM contribution, with a CP-violating NP phase below 10° (see, for 
instance, Ref. [7]). 

The outline of this paper is as follows: in Section [2j we discuss the strategy to extract 
7 and the penguin parameters from the B^ s — > J / ?pKs modes. In Section [3j we have 
a look at the picture for the relevant observables emerging from the current data, and, 
after introducing the framework for the LHCb feasibility study in Section |4| discuss then 
in Section [5] the prospects for the resulting determination of 7. In Section [6| we focus on 
the extraction of the penguin parameters, and discuss then the control of their impact 
on the measurement of (sin2/3)j/^ s in Section[7j Finally, we summarise our conclusions 
and give a brief outlook in Section |8j 



2 The Basic Strategy 



The Bg — > J /ipK$ transition originates from a b — > ccd quark-level process and has a 
decay amplitude of the following structure: 
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where A^ and Ap denote strong amplitudes that are related to the "tree" and "penguin" 
topologies (with internal j G {u, c, t} quarks), respectively, which are illustrated in Fig. [TJ 
On the other hand, the 

= V qd V q * b (2) 

factors contain the relevant combinations of elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. If we make use of the unitarity of the CKM matrix and apply the Wolfen- 
stein parametrisation |8) as generalised in Ref. M, we can write ([!]) as 
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are CP-conserving, "hadronic" parameters, where 
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with A = \V US \ = 0.22521 ± 0.00083 are CKM parameters. It should be noted that R b is 
one of the sides of the UT. 

The Bg — > J/ipKs channel is a decay into a CP eigenstate with eigenvalue —1, and 
offers a time-dependent CP asymmetry of the following structure: 



a C p(t) = 



T(B° S -» J/^jK s ) - T(B° S -» J/xpK s ) 

r( jB o ^ j/^ Ks) + r( J eo j/^ Ks) 
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(7) 

cosh(Ar s t/2) - A A r(B° s ->■ J/^-fsTs) sinh(Ar s t/2) ' v ; 

where C(S° — >■ J/ipKs) and S^-Bj? — > J/ipKs) describe direct and mixing-induced CP 
violation, respectively, and A&r(Bs — > J/ipKs) is accessible thanks to the width differ- 
ence Ar s of the 5 s -meson system. In terms of CP-violating phases and the hadronic 
parameters introduced above, these observables read as follows: 
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where <p s is the CP-violating B° s -B® mixing phase, which is in the Standard Model (SM) 
given by <^ M = -2A 2 r/ = -(2.12 ± 0.11)° but may well be enhanced by NP effects Q. 

It should be noted that these quantities are not independent from one another, sat- 
isfying the relation 

[C(B° S -> J/^K S )} 2 + [S(B° S J/^K S )} 2 + [A AT (B° S -> J/^s)] 2 = 1. (11) 

However, we have still another independent observable at our disposal, which is the 
CP-averaged rate (2] 

(T(B° S -> J/^s)} = $ J7 ^ S x |AT| 2 x [1 - 2acos£cos 7 + a 2 ] , (12) 

where 3>j/^ Ks denotes a phase-space factor, and Af = — AA 

In the case of the decay B® —¥ J/ipK s , the counterpart of the penguin parameter 
ae td enters with the doubly Cabibbo-suppressed parameter e = A 2 /(l — A 2 ) = 0.053 as 
follows: 
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where the hadronic parameters are defined in analogy to the — > J/ipKs case [2], 
and the primes remind us that we are dealing with ab-^s mode. The corresponding 
observables can be obtained from the expressions given above by simply making the 
following replacements: 

a^-a', 9^9', A ->■ A', N W = - ^-Af. (14) 

V e A 

Concerning the mixing-induced CP asymmetry S(B® — > J/ipKs), we have to replace, in 
addition, the B° s -B° s mixing phase by its 5^-meson counterpart 0^, which is in the SM 
given by <^ M = 2(3. 



Thanks to the e suppression of the penguin parameter in (13), the B® — > J/ipKs 
decay is referred to as the "golden" channel for the measurement of <j)d- In the 5-factory 
era of the last decade, it was justified in terms of accuracy to simply neglect this term. 
On the other hand, we are entering a new territory for precision B physics in the LHC 
era which makes it mandatory to take the penguin effects into account. This was already 
anticipated in 1999 in Ref. (2], where it was proposed to use the B® — > J/ipKg channel 
to accomplish this task and to extract the angle 7 of the UT, exploiting the feature 
that the decays B° s ->- J/ipK$ and B Q d ->- J/ipKs are related to each other through an 
interchange of all down and strange quarks, i.e. by the [/-spin flavour symmetry of strong 
interactions. The relevant observables are the CP asymmetries C{B° S — > J/ipKs) and 
S(B® —7- J/ipKs) as well as the following ratio of the CP-averaged rates: 
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The [/-spin flavour symmetry implies 



9 = 9\ (16) 



as well as 

A' = A. (17) 

The three observables depend then on 7 as well as a and 9. Consequently, we can convert 
the measured values of C(B® — > J/ipKs), S(B® — > J/ipKs) and H into 7, a and 9. In 
Sections [5] and |6j we discuss in detail how the implementation of this strategy works in 
practise at the LHCb experiment. 

3 Picture from Current Data 

The CDF collaboration has recently reported the observation of the B® — > J/ipK^ chan- 
nel, with the following result [I]: 

BR(E° -> J/ibK 3 ) , , ,„ . , N 

BR(E° J/^K) = ± °- 0070 ( stat -) ± 0-0041(syst.) ± 0.0050(frag.), (18) 
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where the last error is associated with the ratio f s /fd of the B s and Bd fragmentation 
functions. Using the value BR(5° -> J/ipK°) = (8.71 ±0.32) x 10" 3 of the Particle Data 
Group (PDG) (10| yields 



BR(B° S J/ipK ) = (3.53±0.61(stat.)±0.35(syst.)±0.43(frag.)±0.13(PDG)) x 10~ 5 . 

(19) 

It is interesting to compare this result with the branching ratio of the B® — > J/i/jtt° 
channel. If we neglect penguin annihilation and exchange topologies, which contribute 
to B® — > J/ipn but have no counterpart in B° s — > J/ipK° and are expected to play a 
minor role (which can be probed through B° s — > J/ipn ), the SU(3) flavour symmetry 
applied to the spectator d and s quarks implies 



BR(£?° -» J/j,K°) r Bd su(3) 
2BR(B^J/^)r Bs ^ J/ipKs 



(20) 



where the factor of 2 is associated with the wave function of the 7r°, whereas the T Bi 
and tb s denote the Bd and B s lifetimes, respectively. Using the PDG values T~B d / tb s = 
1.036 ± 0.018 and BR(B° d J/^vr ) = (1.76 ± 0.16) x 10~ 5 , we obtain 



—ST/(3) 



1.01 ±0.25, 



(21) 



which agrees very well with the SU(3) expectation in (20), although the error is still 
sizable. 

In the following analysis, we will use the hadronic parameters extracted from the CP 
violation data for the B® — > J/ipir channel as a guideline, using the numerical values 
obtained in Ref. |3|. Making the same assumptions as in ( J20| , we have then 

(22) 



a G [0.15, 0.67], 9 G [174, 213]°, 



which correspond to 7 = (65 ± 10)°. The numerical value in (21) gives us confidence in 
this procedure. The observables of the B® — > J/ipKg channel then read as 



and 



H = 1.53 ±0.31, 



C = C(B° S J/ipK s ) = C(B° d -> J/t/j7i c 



-0.10 ±0.13. 



(23) 



(24) 



Concerning the mixing-induced CP violation S = S(B® — > J/ipKs), we can calculate 
this observable with the help of Eq. (|9| as a function of the B®-B° s mixing phase 4> s for 
given values of H, C and 7. The corresponding formulae are given as follows: 
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with 



_ 2 [(1 + cH) sin 7 cos 7 ] 2 - (1 - H)(l - e 2 H) sin 2 7 - e [(1 + e)HC cos 7 ] 2 
~ [(1 - e 2 H) sin 7] 2 + [e(l + e)HC cos -f} 2 

= [(1 + e)ff C cos 7 f + [(!-#) sin 7 f 
[(1 - e 2 #) sin 7 ] 2 + [e(l + e)FC cos 7 ] 2 ' 



In these expressions, we have used the [/-spin relation (16) but did not make any other 



approximation. If we assume the SM, we arrive at the following prediction: 

S(B° J/^ s )|sm = 0.54182*1* igSloia = 0.54183. (30) 

As we have noted above, NP may well enter the B® — > J/ipKo, channel through 
contributions to B° s -B° s mixing, yielding a value of 4> s different from 0f M . The most recent 
compilation of the CDF and D0 collaborations using measurements of CP violation 



in Bg —7- J/ip(f> [IT] can be found in Refs. 12 and 13 , respectively. Unfortunately, 
the situation is not conclusive, although the CDF and D0 data are consistent with 
each other: the CDF collaboration finds the ranges <p s e [-59.6°, -2.29°] ~ -30° V 

[-177.6°, -123.8°] 150° (68% C.L.), while the D0 collaboration gives a best fit 

value around <p s ~ —45°, taking also information from the dimuon charge asymmetry 
and the measured B® — > Di*^ + D^ branching ratio into account. 

Interestingly, we can also probe this kind of NP through the mixing-induced CP 
violation in B® — > J ft/) Kg. In Fig. [2J we show the correlation between S(B® — > J/ipKs) 
and sin0 s , which can be determined from the time-dependent angular analysis of the 



Bg — >• J/i/j(f) channel 11 . In this figure, we assume that NP is only present in B° s -B° s 
mixing and have used the formulae give above to calculate curves for different values of 
H and 7. We observe the following features: 

• Thanks to the penguin contributions, we can distinguish between the twofold am- 
biguity for cj) s originating from the measured value of sin0 s , in particular also 
between <p s ~ 0° (as in the SM) and <p s ~ 180°. 

• For cf) s ~ —30° and (f) s ~ —45°, the mixing-induced CP asymmetry is ~ while it 
approaches —1 for <p s ~ —150°. Consequently, we could then also clearly identify 
the scenarios corresponding to the CDF and D0 data discussed above. 

In view of these observations, already first experimental information on the mixing- 
induced CP violation in B® — > J/ipK^ would be very interesting to complement the 
analyses of CP violation in B® — > J/ip<p. Needless to note that this feature relies on 
sizable penguin contributions to the B° s — > J/ipK^ channel. 



4 Framework for the LHCb Feasibility Study 

With an anticipated integrated luminosity of 1 fb _1 by the end of 2011, the LHCb 
experiment expects to be able to measure sin0 s with a precision between 0.04 and 
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Figure 2: Correlation between sin(j) s and 5(5° — > J/ipK$) for different values of H (left 
panel, with 7 = 65°) and different values of 7 (right panel, with H = 1.53). The impact 
of a variation of C = —0.10 ± 0.13 is small. 



0.07 [14]. The B® — > J/ipK§ decay should then already be clearly visible with the 
following uncertainties: 

BR(5° -> J/ipKs) , if, \ 

B R(Bg ->• J/^Ks) = °'° 403 ± °- 0065 ( stat -) ± 0-0029(frag.), (31) 

based on an estimated event yield of 60 000 B\ -»■ J/^K S and 700 B° s -»■ J/^if s decays, 
and where the error on the fragmentation fractions is based on the method proposed 



in Ref. 15 . Concerning the CP-violating observables, we expect a sensitivity of AC = 
AS = 0.34, which might allow us to get first information on these quantities. 

In order to explore the feasibility of the B° s — > J / ipKs measurement at LHCb with 
higher integrated luminosities, we perform a toy Monte Carlo (MC) analysis. For this 
study we assume the following situation: Already with an integrated luminosity of lfb -1 
we can pin down the B° s -B° s mixing phase. Hence, for the two milestone scenarios it will 
be well known. In the following discussion we will consider the SM case (j) s = —2.12°. 
Secondly, the used signal and background yields are extrapolated from Ref. (HI. The 
number of B® — > J/ipK s events listed there is scaled, using the measured B® — >■ J/ipK s 
branching fraction, to give an expected yield of 4000 and 70 000 B° s — > J / tpK$ events for 
the integrated luminosities of 6 fb _1 and 100 fb _1 , respectively. These event yields are 
based on an estimated inclusive J/ ip cross section of cri nc ijM = 262 fib at y/s = 14 TeV, 
the LEP fragmentation functions published by HFAG (l6|, and a bb cross section at 
yjs = 14 TeV of a b i = 568 fib, which was extrapolated from the recently measured bb 



cross section of a b i = (284 ± 20 (stat.) ± 49 (syst.)) fib at y^s = 7 TeV 17 . 

For the MC analysis the reconstructed mass and lifetime distribution of both signal 
and background are modelled. The former is used to obtain the errors on the CP 
observables C and S, while the latter is used for the determination of the statistical 
error on H. The error on H coming from the ratio f s /fd of fragmentation functions, 
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Parameter 


Value 


Parameter 


Value 




568 fib 


B s lifetime resolution 


0.041 ps 




262 fib 


mistag fraction 


0.36 


fd/ fs 


3.55 ±0.26 


B mass resolution 


15.2 MeV/c 2 



Table 1: Overview of the parameters used by the toy MC. 



affecting H via (18), is again determined using the strategy proposed in Ref. 15 . The 
most import parameters used by the toy MC to describe both distributions are listed in 
Table E 

The modelled B s lifetime distribution, given in the left panels of Fig. |3j has three ma- 
jor components: the B® — > J/i])K§ signal (plotted for B s and B s combined), an inclusive 
J/ijj background and a combinatoric bb background. The B s and B s signals are described 
individually using the time evolution equations (11) and (12) in Ref. [2], respectively, 
and include a dilution factor due to mistagging of the original B meson. The inclusive 
J/ijj sample contains both prompt and non-prompt (e.g. from B decays) J/ip candidates, 
which results in an asymmetric distribution. The combinatoric bb background describes 
all other background from bb events; other sources of background are considered negli- 
gible. All lifetime distributions are convoluted with a Gaussian resolution model, which 
causes some events to have a negative lifetime. The resolution for the B° s — > J/ipKs 
signal is a = 0.041 ps jH]. 

The modelled mass distribution, also given in Fig. [3j describes the two B mass peaks 
(Bd on the left, and B s on the right and in the insert) and an exponential background 
for both the inclusive J/ip and combinatoric bb sample. 

The signal yield is small compared to the number of background events, as can be 
seen from the inserts in Fig. |3j making it hard to distinguish the B s mass peak from 
the background. This is the main reason for choosing the high luminosity scenarios of 
6 fb _1 and 100 fb _1 . Nonetheless, the lifetime distribution shows that the background 
is predominantly prompt. The long lifetime region (t > 2 ps) is therefore sufficiently 
sensitive to B° s — > J/i/jKg in order to extract C and S. This is confirmed by the plot 
of the time-dependent CP asymmetry acp(t) defined in equation ^ and given in Fig. [5] 
for our two milestone scenarios. 

The toy MC study is mainly used to estimate the errors on the relevant observables. 
The actual values of these observables are determined by their dependency on a, 9 and 
7, given in equations (JsJ) , ^ and (15). The input values for a, 9 and 7 were taken from 
the B^ — > J/ipir analysis whose result was already given in (22). We take the central 
values of these intervals, 

a = 0.41 , 9 = 194° , 7 = 65°, (32) 

which result in the following picture for the statistical uncertainties at LHCb: 

6 flT 1 : AH = 0.069, AS = 0.14, AC = 0.14, (33) 

100 ftT 1 : AH = 0.052, AS = 0.035, AC = 0.035. (34) 

In principle these statistical errors could depend on the input value for a, 9 and 7, but 
we will assume them constant for the remainder of this paper. 
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Figure 3: The lifetime distribution for an originally tagged B s meson (left) and the total 
mass distribution (right) for the decay B° s — > J '/if)Ks and its background, as used in the 
MC analysis. This shows good sensitivity to the signal in the long lifetime region. The 
number of events correspond to a data sample of 6 fb -1 (top) and 100 fb -1 (bottom) 
respectively. 
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Figure 4: Toy MC simulation of the time-dependent CP asymmetry acp{t) denned in 
equation (J7| for the 6 fb _1 (top) and 100 ffcr 1 (bottom) scenarios. 
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For a more detailed discussion of the experimental feasibility study, the reader is 
referred to Ref. [l~8| . 



5 Extraction of 7 at LHCb 

First information about 7 could in principle be obtained from H as the general structure 



of this observable in (15) implies the following inequality 19 

H > [l -2ecos 2 7 + C(e 



sin 2 7 



(35) 



However, as the experimental data indicate a value of H > 1, the corresponding bound 
on 7 is not effective. 

As discussed in Ref. j2], using the B®-B® mixing phase as an input, we can convert 
the B° s — > J/i/jKs CP asymmetries C and S into a contour in the 7-a plane, which 
is theoretically clean. Using the t/-spin symmetry, we can determine another contour 



from H and S(B® — > J/if)K$). Thanks to the [/-spin relation in (16), the intersection of 
these curves allows us then to extract 7 and the penguin parameter a. Knowing these 
parameters, the CP-conserving strong phase 9 can be extracted as well. In Fig. |5j we 
illustrate the corresponding situation in the 7-a plane for our LHCb study. 
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Figure 5: Determination of 7 and a from intersecting contours. Situation resulting from 
our LHCb feasibility study for 6 fb _1 (left) and 100 fb _1 (right). 

The angle 7 can also be extracted, together with a and 9, from a simultaneous x 2 fit 
to C, S and H using the expressions in (|8]) , ([9]) and (15). The corresponding numerical 
results read as follows: 



100 ftT 1 : 



: 7 = (65 ±10)°, 
7= (65.0 ±3.2)°, 



= 0.410 ±0.079, 
0.410 ±0.023, 



9 = (194 ±17)°, (36) 
9 = (194.0 ±4.1)°. (37) 



Here we give only the statistical errors, which show that LHCb should have sensitivity to 
extract 7 from the [/-spin-related B s ^ — > J/ipKs decays. Using other avenues to extract 
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7, in particular theoretically clean strategies employing B® — > D Z ^K ± and similar modes, 
LHCb expects the following accuracies for 7: 

6 ftr 1 : 7 = (65.0 ±3.1)° [14], 100 fb" 1 : 7 = (65.0 ± 1.0)° §, (38) 

which should be compared with the current situation emerging from B -> D^K^ 
analyses: 

7 = (71ig)° (CKMfitter M), (78 ± 12)° (UTfit Ell). (39) 



Consequently, LHCb should make significant progress with the determination of 7 with 
respect to the unsatisfactory current situation. 

Concerning the B s ^ — > J/ipKs strategy, we find a good sensitivity for this angle in 
our feasibility study. We have also to deal with theoretical uncertainties originating from 



[/-spin-breaking effects. As far as the corrections to (16) are concerned, they actually 
play a very minor role for the determination 7 as the a' terms in (15) are strongly 
suppressed by the tiny e parameter in [2]; numerically, they give a correction for 7 well 
below 1°. However, ?7-spin-breaking effects enter also through the ratio |*4'/*4|, which 
affects the determination of H as defined in ((TBI. We find that an uncertainty of 10% 



for the relevant form-factor ratio F B o K o(Mj / ^;l ) /F B o^o (Mj /ib ; 1 ), as considered in 



Ref. 13], corresponds to an error of about A7|ff = 12°, which essentially matches the 



statistical error of 10° for our 6 fb _1 scenario in (36 ). For an LHCb upgrade with 100 fb -1 , 
a form-factor uncertainty of 2.5% would be required in order to match (37), which serves 
as the ultimate benchmark for future calculations of this quantity. We hope that already 
by the end of the second LHCb data taking period, i.e. 2014-15, we will have a better 
understanding of 577(3) breaking through studies of various other [/-spin- and S77(3)- 
related I?-meson decays, as well as through improved theoretical studies, including in 
particular also non-factorisable [/-spin-breaking corrections to |.A'/*4|. The current data 
from the Tevatron do not indicate sizable non-factorisable S77(3)-breaking effects |22|, 



which is also supported by (21). 



An important aspect of the determination of 7 through the B s ^ — > J/i/iKq strategy 
is that it may reveal a NP contribution to the B° s — > J/ijiKg decay amplitude. This 
would be indicated by a value of 7 in conflict with the clean determinations through 



the pure "tree" decays summarised in (38). In the following discussion, we assume that 
we have negligible NP effects of this kind, i.e. assume the SM expressions for the decay 
amplitudes in ^ and (13). 



6 Clean Determination of the B® — > J/iJjKs Penguin 
Parameters at LHCb 

The major application of the B® — > J/ipKs channel at LHCb will be the extraction of the 
hadronic penguin parameters (a, 6) and to take them into account in the determination 
of the B®-B® mixing phase <p d from (sin 2(3)j/^ Ks . Since our goal is to minimise here the 
[/-spin-breaking corrections, we shall refrain from using the H observable and will use 7 



as given in (38) as an input. Following Ref. 3 , we can then convert the CP asymmetries 



S and C into theoretically clean contours in the 9-a plane; their intersection allows us 
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Figure 6: Determination of a and 9 from intersecting contours. Situation resulting from 
our LHCb feasibility study for 6 fb" 1 (left) and 100 fb" 1 (right). 



then to extract a and 9. In Fig. [6j we illustrate the corresponding situation arising 
from our LHCb feasibility study. In order to guide the eye, we show also the contour 
corresponding to H, which would be interesting to obtain insights into ^/-spin-breaking 
effects. Performing again a \ 2 fit yields the following results: 

6 ftT 1 : a = 0.41 ± 0.14 (stat.) ± 0.02 (7) , 9 = [194 ± 16 (stat.) ± 0.8 (7)] , (40) 
100 ftT 1 : a = 0.41 ± 0.03 (stat.) ± 0.01 (7) , 9 = [194.0 ± 4.0 (stat.) ± 0.3 (7)] .(41) 

We observe that LHCb is expected to be able to determine the B° s — > J/ipKs penguin 
parameters with impressive precision. It should be emphasised that these parameters 
rely only on the SM expression for the decay amplitude and are theoretically clean. 



7 Controlling Penguin Effects in the Measurement 
of (sin 2/3) j/ipKs a ^ LHCb 

The values of the penguin parameters a and 9 are interesting from the point of view of 
strong interaction dynamics. However, they have also an important practical applica- 
tion for testing the Kobayashi-Maskawa mechanism of CP violation. Using the [/-spin 



relation (16), we can relate the B® — > J/i/jKs parameters to their B® — > J/ipKs counter- 
parts, which allows us then to control the hadronic uncertainty of the measurement of 
the B®-B® mixing phase (pa through the mixing-induced CP violation in B® — > J/ipK^. 



The corresponding, generalised expression reads as follows [3 

S(B° d -». J/^ s ) 



v/1 - C(B° d J/Ws) 



sin^ + A^), (42) 
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Figure 7: Determining the penguin correction Atfid to sin (2/3) from intersecting contours. 
Situation for 6 fb _1 (left) and 100 fb" 1 (right). 



where 



sin A(p d 
cos Aifid 



2ea cos 9 sin 7 + e a sin 27 



Ay 1 - c(j/^K s y 

1 + 2ea cos cos 7 + e 2 a 2 cos 27 



(43) 
(44) 



with N = 1 + 2ea cos 9 cos 7 + e a , so that 



tan A0d 



2ea cos 6 1 sin 7 + e a sin 27 



1 + 2ea cos 9 cos 7 + e 2 a 2 cos 27 
The current value of the mixing-induced CP violation is given by 

S(B° d -»■ J/V>K°) = 0.655 ± 0.024, 

while the direct CP asymmetry is currently measured as 

C(B° d J/1PK ) = -0.003 ± 0.020. 



(45) 



(46) 



(47) 



These numbers are averages of the BaBar and Belle data over the B® — > J/ipKg and 
B d — > J / ipKi, final states, as obtained by the Heavy Flavour Averaging Group [16] . The 



value in ( [47] ) implies that the deviation of the square root in (42) from one is at most at 
the level of 0.0002, i.e. completely negligible. 

Using these formulae, we can convert the contour plots in Fig. [6] into curves in the 
9-A(j)d plane, as shown in Fig. [7] In this analysis, we have uncertainties from [/-spin- 
breaking corrections to (16), which we parametrise as follows: 



to! 



9 = 9' + A9. 



(48) 
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Figure 8: Correlation between the hadronic shift A0^ and its error for our LHCb feasi- 
bility study for different input values of the penguin parameters a and 9. Situation for 
6 fb _1 (left) and 100 fb _1 (right); note the different scales of the y axes. 



If we vary f G [0.85, 1.15] and A6> G [-20°, 20°], and perform a x 2 fit to C and S, we 
obtain 

6 fir 1 : A<p d = [-2.23 ± 0.78 (stat.) ± 0.12 ( 7 ) ±° f 7 (A6)] ° , (49) 

100 fb _1 : A<j) d = [-2.23 ± 0.19 (stat.) ± 0.04 (7) ^(f) t f 7 {A6)}° . (50) 

In order to study the dependence on the input values of the penguin parameters for 
our analysis, we show the correlation between the hadronic shift A0^ and its statistical 
error at LHCb in Fig. |8j The corresponding curves show nicely that we can actually 
well determine the A<pd correction for a wide range of (a, 6) that should contain the 
"true" values of these parameters. It should also be noted that already a small penguin 
contribution with a = 0.1 gives a correction of ~ —0.5°. 

In order to fully appreciate the results of this study, we have to look at the precision 
for (sin2/3) JAWCs at LHCb: we expect 0.022 for 21b" 1 |], and 0.003-0.010 (depending on 
to be determined systematic errors) for 100fb _1 |6j. We extrapolate these numbers into 



0.014 for 6 fb _1 , which should be compared with (46). For a central value of S(B® — > 
J/if)K°) ~ 0.655, these numbers correspond to errors of about 1° and (0.2-0.8)° for the 
phase itself. In order to match these very impressive precisions, we definitely have to 
control the doubly Cabibbo-suppressed penguin contributions in — > J/ipK s . Looking 



at the numerical results in (49) and (50), we observe that we can actually achieve this 
goal. Already with 6fb~ 1 we have to control the penguin effects, for an upgrade of LHCb 
this is absolutely necessary to fully exploit the tremendous statistics. Measurements 
along these lines may eventually allow us to resolve CP-violating NP contributions to 
B^-Bj mixing. 

Let us finally note that we can - in addition to the direct CP violation in the 
Bj — > J/ipKs t L channels |3| - add another observable to this analysis. It is the direct 



15 



CP asymmetry of the — » J/ip7r + decay, which we get from the B° s — > J/ipKg topolo- 
gies in Fig. [T] by replacing the strange spectator quark by an up quark. The resulting 
contours in Fig. [7] are similar to those related to C. Already with 1 fb _1 , as expected 
by the end of 2011, we find from a feasibility study similar to that described above that 
the .B+ — > J/ip7T + mode should be clearly visible at LHCb, with a statistical error of 
C(B^ J/ipiT^) = 0.005. However, the production asymmetry between B and B states 
at the LHC (which, unlike the Tevatron, is a pp collider) is expected to be non-negligible 
and will have to be measured and corrected for. This is likely to induce irreducible sys- 
tematic errors at a few per-cent level for the measurement of C(B^ — > J/-07T ), which 
have to be studied in further detail. 

8 Conclusions 

As was pointed out about one decade ago, the decay B° s — > J/ipKs is the £/-spin partner 
of the "golden" mode B® — > J/ipKs and allows an extraction of the UT angle 7 and of 
the corresponding penguin parameters. In the summer of 2010, the CDF collaboration 
has announced the first observation of the B Q S — > J/ipRg transition, with a branching 
ratio in agreement with an SU(3) relation to the B® — > J/ipn decay. 

We have calculated a correlation between the mixing-induced CP violation in the 
B° s — > J/if)Ks channel, S, and sin0 s , which serves as a "target region" for the first CP 
violation measurement in this decay. We observe an interesting pattern for values of the 
B° s -B° s mixing phase (f) s in the region of the current measurements from the Tevatron: 
we can distinguish between the twofold solutions for <ft s arising from the measurement of 
sin 0s, with S ~ and S ~ —0.8 for (fi s ~ —30° and (f> s ~ —150°, respectively, whereas 
we expect S ~ +0.5 in the SM. Thanks to improved measurements of input parameters, 
these predictions can be sharpened in the future. 

By the end of 2011, i.e. the end of the first LHCb data taking period, we expect a 
clear signal for B° s — > J/ipK$, with very crude first information on the corresponding CP 
asymmetries. Extrapolating from published LHCb studies, we have performed a detailed 
feasibility study for the measurement of the B° s — > J/ipKs observables, yielding about 
4000 and 70000 signal events for 6 fb _1 and 100 fb _1 , respectively. These integrated 
luminosities refer to the end of the second LHCb data taking period of 2014-15 and a 
subsequent LHCb upgrade. Using input parameters as suggested by current data, we 
have shown that the B s d — > J / ipKo, strategy offers indeed another determination of 7 for 
LHCb. However, we conclude that the most important application of this method will be 
the determination of the corresponding hadronic penguin parameters and their control 
in the measurement of (sin 2/3)j/^k s , which will allow us to match the corresponding 
experimental precision. Studies along these lines may eventually allow us to resolve NP 
in Bj-Bj mixing. 

Acknowledgements 

We would like to thank Marcel Merk and Wouter Hulsbergen for carefully reading the 
manuscript and useful comments and suggestions. 



16 



References 

[1] CDF Collaboration, CDF Note 10240, August 2010, see also 



http : / / www-cdf.fnal.gov / physics / new /bottom /100708 .blessed-Bs JpsiK / 



[2] 
[3] 

[4] 
[5] 
[6] 
[7] 
[8 
[9 

[10 
[11 

[12 

[13 

[14 
[15 

[16 

[17 
[18 
[19 



R. Fleischer, Eur. Phys. J. C 10 (1999) 299 |arXiv:hep-ph/9903455|. 

S. Faller, M. Jun g, R. Fleischer and T. Mannel, Phys. Rev. D 79 (2009) 014030 
|arXiv:0809.0~842l [hep-ph]]. 

M. Ciuchini, M. Pierini and L. Silvestrini, Phys. Rev. Lett. 95 (2005) 221804. 

LHCb Collaboration, Report CERN-LHCC-2003-030 (2003). 

LHCb Collaboration, Expression of Interest CERN-LHCB-2008-019 (2008). 

P. Ball and R. Fleischer, Eur. Phys. J. C 48 (2006) 413 |arX iv:hep-ph/0604249| . 

L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945. 

A. J. Buras, M. E. Lautenbacher and G. Ostermaier, Phys. Rev. D 50 (1994) 3433 
|arXiv:hep-ph/9403384| . 

K. Nakamura et al. (Particle Data Group), J. Phys. G 37 (2010) 075021. 

I. Dunietz, R. Fleischer and U. Nierste, Phys. Rev. D 63 (2001) 114015 |arXiv:hep-| 
|ph/0012219]; A. S. Digh e, I. Dunietz and R. Fleischer, Eur. Phys. J. C 6 (1999) 647 
|arXiv:hep-ph/9804253| . 

CDF Collaboration, CDF/ANAL/BOTTOM/PUBLIC/10206, July 2010, see also 



http : / / www-cdf . fnal .gov / physics / new /bottom 



D0 Collaboration, D0 Note 6093-CONF, July 2010, see also 



http : // www-dO . fnal .gov / Run2P hysics /WWW / results /b . htm 



LHCb Collaboration, |arXiv:0912.4179 [hep-ph]] 



R. Fleischer, N. Serra and N. Tuning, Phys. Rev. D 82 (2010) 034038 
|arXiv: 1004.3982] [hep-ph]]. 

E. Barberio et al. (Heavy Flavour Averaging Group), |arXiv:0704.3575 [hep-ex]]; for 



the most recent updates, see http://www.slac.stanford.edu/xorg/hfag, 



R. Aaij et al. (LHCb Collaboration), |arXiv:1009.2731| [hep-ex]. 
K. De Bruyn, Master thesis, University of Amsterdam, in progress. 



R. Flei scher and S. Recksiegel, Phys. Rev. D 71 (2005) 051501 larXiv:hep- 
| ph/0409137| ; R. Fleischer and T. Mannel, Phys. Rev. D 57 (1998) 2752 |arXiv:hep- 
ph/9704423|. 



17 



[20] CKMfitter Collaboration, |http://ckmfitter.in2p3.ir7 



[21] UTfit Collaboration, http://www.utfit.org/UTfit/ 



[22] R. Fleischer, Eur. Phys. J. C 52 (2007) 267 |arXiv:0705.1121 [hep-ph]]. 



18 



